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ABSTRACT 

Context. The study of the evolution of the Tully-Fisher relation has been controversial in the past years. The main difficulty is 
in determining the required parameters for intermediate and high redshift galaxies, given the cosmological corrections and biases 
involved. 

Aims. This work aims to identify the main problems of the study of the Tully-Fisher relation at high redshift using optical emission 
lines, in order to draw conclusions about the possible evolution of this relation in the B, R, and I-bands. 

Methods. With this aim, the rotational velocities obtained from the widths of different optical lines using DEEP2 spectra are compared. 
This method allows reaching higher redshifts against the rotation curve one, due to spatial resolution limitations. Morphology has been 
determined via HST images, using and comparing different methodologies. Instrumental magnitudes are then corrected for K and 
extinction and the absolute magnitudes derived for the concordance cosmological model. Finally, the optical Tully-Fisher relations in 
B, R, and I-bands at different redshifts up to z = 1.3 are derived. 

Results. Although most studies (this one included) find evidence of evolution, the results are not conclusive enough, since the possible 
luminosity evolution is within the scattering of the relation, and the evolution in slope is difficult to determine because at high redshift 
only the brightest galaxies can be measured. Nevertheless, our study shows a clear tendency, which is the same for all bands studied, 
that favours a luminosity evolution where galaxies were brighter in the past for the same rotation velocity. This result also implies 
that the colour of the Tully-Fisher relation does not change with redshift, supporting the collapse model versus the accretion model 
of disc galaxy formation. 
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1. Introduction 

The Tully-Fisher relation (hereafter TFR) relates the lumi- 
nosity to the maximu m rotational velocity of spiral galaxies 
dTullv & Fishet] 1 1977b. This relation is va lid for all types of 
spirals. For example, IZwaan et"al I (119951) obtained the same 
TFR for high and low surface-brightness galaxies, while 
ICourteau et alj d2003l) came to the same conclusion for barred 
spirals. 

Then, the TFR represents a relation between fundamental 
galaxy parameters, such as its total mass, and the mass locked in 
stars. This valuable info rmation about the relation between dark 
and luminosity matter (IZwaan et al.ll 19951) has important impli- 
cations for the study of the evolution of disc galaxies. Its slope is 
an input parameter for galaxy formation models, and its absolute 
luminosity level provides constrai nts for the mass-to-light ra- 
tios a nd the initial mass functions dPortinari & Sommer-Larsenl 
120071) . Finally, the TFR is an important distance estimator that 
has bee n traditionally used f or measuring the Hubble constant 
H (e.g. lTullv & Pierc3l2000l) . 

For these reasons, the study of the possible evolution of the 
TFR by comparing what has been obtained from galaxies at high 
redshift with their local counterpart, has outstanding importance 
and far reaching consequences for determining fundamental cos- 
mological parameters in the st udy of structure formation, and in 
the evolution of disc galaxies dNavarro & Steinm etz 2000). 

In the past years, some groups have measured the TFR at 
different redshifts. Pioneer studies, most of them in the B- 



band, seem to show evolution, although there was no consen- 
sus on wheth e r this evolution was on slope, zero point, or both. 
Zieg ler et al.1 (120021) find a shallower relation at high redshift 
(z=l) than what has been measured in local samples. In their 
study, they claim that less massive galaxies were 1 or 2 magni- 
tudes brighter in the past (at fixed rotational velocity), while the 
most massive would follow the local relation. Bohm et al. (2004) 
find the same evolution on slope and evidence for a luminosity 
evolution with a look-back time of AMb ~ -1 magnitudes at 
redshi ft z=l. However, other groups, such a s Rix et al. [jl997h 
or lBamford. Aragon-Salamanca & Milvang-Jensenl (12006). find 
luminosity evolution but no slope change. Moreover, these au- 
thors find differences in luminosity evolutions ranging from 0.2 
to 2 magnitudes. 

Since optical bands, especially in the blue, are more affected 
by extinction, infrared TFR have been studied since very early 
on (lAaronson. Huchra & Mouldll 1979b . When t he infrared TFR 
was re cent ly applied to evol utionary studies, IConselice et al.l 
(120051) and lFloresetail d2006l) did not find a K-band T FR evo- 
lution. However, in a recent work, iPuech et ail (120081) do find 
evolution in the sense that galaxies had been fainter in the past, 
a resul t opposite to what is found in the optical bands (for ex- 
ample, Bohm et al. 2004). These results suggest that the TFR 
could be dependent on the observed photometric band. An al- 
ternative method for studying the TFR has been derived by 
iBell & De Jongl d200lh using the stellar-mass obtained from the 
mass-to-light ratio. These authors claim that it has the advan- 
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tage of being band-independent. The principal disadvantages 
are, however, its larger scatter, and an initial mass function (IMF) 
must be assumed. Even so, some studies have used this method 
(e.g.. lConselice et aT]|2005l) . finding no evolution in the stellar- 
mass TFR. Nevertheless, from t he combination of t he rotational 
velocity and velocity dispersion, Kassi n et al. I d2007h reduced the 
stellar mass TFR scatter, concluding that there is no evolution up 
to z= 1 .2. In summary, all these previous works indicate that the 
issue of TFR evolution if far from being settled. 

While there seems to be a consensus on the absence of 
evolution in the stellar-mass TFR, the results obtained from 
intermediate-redshift optical studies (0.3<z<0.7) show discrep- 
ancies that could be attributed to observational biases. Some 
works have attempte d to pinpoint these limitations, such as the 
one by iBarton et al.l d200ll) . who conclude that the possible ex- 
planations for the TFR outliers (i.e. objects with substantial de- 
viations from the local TFR) could be gas infall, non-uniform 
or trun cated emission, and disc-en hanced star formation. In ad- 
dition, |Kann^pani&BOTto^ d2004l) find that the local TFR out- 
liers associated with kinematic anomalies occupy the same re- 
gion of TFR parameter space as the galaxies respon sible for TFR 
slope evolution found in some high-z studies (e.g. lZiegler et alj 
120021) . which tend to us e mo rphology-blind selecti on criteria. 
Also. lBohm et alj (12004). and lNakamura et all d2006l) conclude 
that the evolution observed is strongl y dependen t on selection 
criteria. Moreover, in a recent study, Vogt (2007) does not find 
evolution up to z — 1.1, when discarding ellipticals and face-on 
spirals from her study, instead of applying the more commonly 
used criteria of selecting galaxies of intense lines or by their in- 
clination angles. Therefore, object selection is of outstanding im- 
portance for drawing conclusive results on the evolution of the 
TFR. However, when extending the study of the TFR to higher 
z, it is difficult to segregate interactive and anomalous galaxies, 
which are frequently included in the sample. 

The rotation velocity measurement is another challenge in 
TFR studies. At high z, the spatial resolution required to resolve 
the rotation curve is comparable to the atmospheric seeing. For 
example, at z — 1 the Galaxy would have an apparent diameter of 
about 3 arcseconds. In these cases, integrated velocities are the 
most convenient way to estimate the velocity field of disc galax- 
ies. From the relation between the velocity field obtained from 
line widths and from rotation curves, Vogt (20071) concludes that 
both methods are equivalent. 

There are other types of errors, such as inclination an- 
gle determination or extinction correction, that tend to shift a 
galaxy along the TFR rather than perpendicular to it dVogtl2007t 
Van Starkenburg et al. 2006); as a consequence, they do not in- 
troduce fake evolution in the TFR at intermediate or high red- 
shift. 

In this work, a study of the rotation velocity from optical 
emission lines has been carried out using different lines, with 
the aim of reliably increasing the redshift range, reaching up to 
Z = 1 .6 in the optical and even more in the infrared. We find that, 
although all optical lines can be used to calculate the rotation 
velocity, the line fitting technique and the signal-to-noise ratio 
play important roles in its proper determination at high redshift. 
In addition, we used three different techniques for morphological 
classification to look for the best method for discarding E/S0 
type and interactive galaxies from our sample. For the first time, 
B, R, and I band TFRs are estimated using a significant sample 
of galaxies up to z = 1.3. With the study of the TFR evolution 
in several optical bands, its possible colour differences versus 
redshift could provide clues to the mechanism responsible for 
this evolution. 



This paper is organised as follows. In Sect. 2, the DEEP2 
spectroscopy, a description of the morphological analysis, and 
sample selection criteria are provided. The study of the opti- 
cal lines and derivation of the rotation velocity and luminosi- 
ties are explained in Sect. 3. Finally the results are given in 
Sect. 4, whereas the last two sections provide the discussion 
of the results and the conclusions. Throughout this article, the 
concordance cosmology with = 0.7, D. m - 0.3 and Ho = 
70 km s _1 Mpc 1 is assumed. All magnitudes are in the AB zero- 
point system. 



2. Data 

The sample consists of galaxies in a 16'xl6' field in the Groth 
Strip Survey (GSS) sky region, which will be observed by the 
OTELO projec t (OSIRIS Tunable Emission Line Object Survey, 
Cepa et al. 200§j). This project will s urvey emission-line objects 
using the OSIRIS (Cepa et al. 2005) tunable filters in the wave- 
length intervals defined by windows through the OH emission 
line forest. OSIRIS is the optical Day One instrument of the GTC 
10.4m telescope. 

The baseline for spectroscopy target pre-selection were the 
galaxies for which DEEP2 spectra (Data Release 3, DR3) in this 
field were available. In the redshift range 0.1<z<1.3, there are 
1200 galaxies in the DEEP2 DR3 sample with a B-band limiting 
magnitude of 24. 



2.1. The DEEP2 Database 



The DEEP2 project dDavis et alJ I2003L I2007I) is a survey us- 
ing the Keck telescopes to study the distant Universe. The sur- 
vey targeted a total of ~10.000 distant galaxies in the redshift 
range 0<z<L4, us ing the DEIMOS multi-object spectrograph 
(Fab er et ai]|2003l) . The grating used was the 1200 1/mm one, 
covering a spectral range of 6500-9100 A with a dispersion 
of 0.33 A/px, equivalent to a resolution R = A/AA 4000. The 
DR3 constitutes the third public data release for this survey and 
includes redshifts plus B, R, and I band photometry, in addi- 
tion to ID and 2D spectra. The photometric data wer e taken 
with the CFH12K mosaic camera dCuillandre et alj200ll) . which 
is installed on the 3.6-meter Canada-France-Hawaii Telescope 
(CFH T). In the DEEP2 photometric catalogue (Data Release 1, 
DR1) dCoil et al. 20041) all magnitudes were corre cted for red- 
dening using Schlegel, Fink beiner & Davisl d 1998b dust maps. 
Magnitude errors from sky noise and redshift are available too. 

As already pointed out, since resolved rotation curves are 
difficult to obtain at high redshift due to the limited spatial res- 
olution, it has been necessary to use the integrated spectra ob- 
tained by fitti ng line profiles to the ID spectra provided by DR3 
dHorndl 19861) . 

2.2. Morphology 

As stated in Sect. 1, morphological classification plays an impor- 
tant role in TFR determinations. Studies with no morphological 
selection criteria include objects such as blue compact galax- 
ies, emission-line SO, or galaxies with peculiar morphologies, 
together with normal spirals. 

In the DEEP2 catalogue, the objects are already divided in 
galaxies, stars, or AGNs, according to a best-fitting template 
technique. However, a morphological classification of the ob- 
jects thus identified as galaxies is still required. As a prelimi- 
nary procedure and with the aim of segregating E/S0 galaxies 
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Fig. 1. Maximum rotational velocity vs. instrumental magnitude, 
in the redshift interval 0.3<z<0.8. Diamonds are galaxies with 
inclination angle lower than 25°. Squares represent galaxies with 
inclination angles between 25° and 75°. Triangles are galaxies 
with inclination angles between than 75° and 80° (the limit in 
our sample). 



from non-E/SO in our sample, we compared the DEEP2 spec- 
tra with a synthetic, early-type galaxy spectrum. The chosen 
template corresponds to a 5 Gyr-old stellar populati on with a 
stella r metallicity of 0.02 Z H (assuming an IMF from ICh abrier 
120031) . which is part of the iBruzual & Charlotl d2003l) collec- 
tion of galaxy templates used by iTremonti et alj d2003h in the 
analysis of SDSS galaxy spectra. Then, the r.m.s. of the differ- 
ences between the template and every object were calculated to 
determine whether they matched each other or not. The result 
was that about one-third of the spectra could be associated with 
E/S0 galaxies. Nevertheless, when the HST images of a statis- 
tically significant subsample of each group were visually anal- 
ysed, the classification mismatch for the non-E/SO group was 
~3%, whereas the one associated with the group of E/S0 galax- 
ies was significantly worse. This result forced us to proceed to 
a visual classification of every galaxy using HST images, as the 
most reliable procedure for an accurate determination of mor- 
phol ogical types . The H ST images used are part of AEGIS sur- 
vey dDavis et al.l 120071) and were obtained with the ACS cam- 
era. There are 809 galaxies in our field that have a DEEP2 spec- 
trum, and V and I-band HST images. In our morphological clas- 
sification we divided the objects into five groups: elliptical/SO 
(8%), spirals (69%), irregulars (6%), interactive (6%), and un- 
known (1 1%). To classify our visually unknown objects, GIM2D 
(Simard 1998) was used. In addition, to somehow estimate the 
reliability of the GIM2D classification, GIM2D parameters were 
also determined for other types. For the galaxies visually classi- 
fied as spirals, the Sersic index was obtained for a random sub- 
sample of 170 galaxies. An apparent classification mismatch was 
found for 40 objects (~24%). However, 25 of them were clearly 
spirals, even though the Sersic index calculated was greater than 
2.5. The Sersic index was also determined for a random sub- 
sample of galaxies visually classified as ellipticals. In this case, 
~60% had a Sersic index lower than 2.5, although half of them 
were clearly ellipticals. These results show that our visual clas- 
sification has clear advantages over that obtained using GIM2D. 
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Fig. 2. Redshift distribution of all DEEP2 galaxies in our field 
(solid line) and of galaxies selected according to our morphology 
classification and inclination criteria (dashed line). 



2.3. Sample selection 

In some studies of the TFR, the sample is sele cted accord- 
ing to the emission line observed. For example, iBarton et al.l 
d2001l) considered only galaxies with H a emission lines, while 
iFlores et all d2006l) studied only galaxies with [Oll]AA 3727 A 
emission. In the present work, all conspicuous optical emission 
lines will be used to expand the redshift range to the maximum 
possible. 

The inclination angle (i) is an usual s election criterion as 
well. For i nstance, [ Nakamura et al. (2006) study galaxies with 
i> 45°, and lKassin et al.l d2007l) choose spirals with 30° <i<70°. 
In the first case, this restriction can include galaxies with large 
inclination errors (almost edge-on galaxies), and in both cases, 
the limits seem too restrictive. With the aim of avoiding selec- 
tion criteria that can bias the results, morphology and inclination 
selection will be carefully studied. The inclination angle can be 
obtained from the major to minor axes ratio as found in the HST 
catalogue. In order to o btain the errors, we used sextractor 
dBertin & Arnoutsll 19961) in the summed V+I images and com- 
pared the results. We estima ted the mean err or to be ~ 2.5 de- 
grees, in agreement with IFlores et alj d2006l) . In Fig. 1, galax- 
ies with three different inclination ranges have been represented. 
No correlation between magnitude and rotation velocity can be 
observed for inclinations lower than 25°. Also, galaxies almost 
edge-on will be more affected for extinction and the error in in- 
clination angle will be larger. 

In Fig. 1, the galaxies with inclination angle up to 75° seem 
to follow a TFR, and in our sample all galaxies have inclina- 
tion angles lower than 80°. For these reasons, only galaxies with 
inclinations between 25° and 80° will be considered. 

Finally, only normal spiral galaxies follow the TFR. Then, 
a first selection of spirals was done using our visual classifica- 
tion. All galaxies belonging to E/S0 and interactive groups were 
discarded. The irregular galaxies were, nevertheless, included, 
since they cannot be distinguished from spiral galaxies at high 
redshift. Either way, the number of irregular galaxies is not sig- 
nificant, and they affect low and high redshift regimes as well. 
There are a few objects visually classificated as unknown and 
without Sersic index that were eliminated because their mor- 
phology was not known. However, it is very difficult to do a 
reliable classification at high redshift, and then peculiar galax- 
ies with kinematic anomalies could be included in our study. 

Figure 2 shows that the redshift distribution of all DEEP2 
galaxies in our field, and that of our final sample, after morphol- 
ogy and inclination selection, are very similar. 
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Fig. 3. Comparison of the velocity widths from five optical lines. In each panel, the solid line shows the least-square fits with slope 
1. All zero-points are within the errors, a) Velocity obtained from [NII]/i6583A and H a , b) comparison between [OIII]/15007A and 
H a , c) velocity from Hg versus [OIII]/15007A, d) comparison between [OIII]/15007A and [OII]^A3727A velocities. 



3. Deriving Tully-Fisher relations 

The final sample of galaxies with redshifts and B, R, and I mag- 
nitudes, classified as spirals, and with inclinations in the range 
between 25° and 80° is composed of 612 galaxies. The next 
step is determining inclination and redshift-corrected magni- 
tudes and rotation velocities. The derivation of the Tully-Fisher 
relation also requires determining absolute B, R, and I-band 
magnitudes. They have been obtai ned by K-correcting B, R , and 
I-band apparent magnitudes using Blanton & Roweis (2007), by 
applying a concordance cosmology for determining the luminos- 
ity distanc e from the redsh ifts, and by correcting for intrinsic 
extinction dTuUv et alj|1998l) . 

3.1. Rotation velocity estimation 

The local TFR has been historically obtained from radio mea- 
surements, specifically from 21cm line widths. The HI emis- 
sion spreads further than optical emission, so that optical and 
radio measurements can be compared, as long as both extend 
up to the flat region or to the maximum of the rotation curve. 
Moreover, the existence of three types of rotation curves, de- 
pending on the relation between the maximum velocity and 
the velocity of the flat region (IVerheiienll200lh . can compli- 
cate the comparison. Also, turbulence plays a role in the op- 
tical versus radio velocity width det ermination. For example, 
Math ewson. Ford & Buchhornl dl992[) compared the projected 
rotation velocity measured from H a rotation curves with the ve- 



locity measured from integrated HI profiles at the 50% level. 
They obtain a difference of 10 km/s, which is attributed to HI 
widths that measure rotational velocity plus internal galaxy tur- 
bulence. This turbulence is important in the most external region 
of the galaxy where the gravity is lower. The optical emission 
does not spread until this region, so the turbulence is negligible 
against the rotational velocity determined using these lines: 

AAc 

max — ~ : r ■ v*) 

i sin(i) (1 +z) 

To measure rotation velocity, we used equation (1), where i is 
the inclination angle, Aq the line central wavelength at z=0, and 
AA the line width at 20% of peak intensity. We determined AA 
with the Gaussian fi tting routines in the Starlink package dipso 
(iHowarthet all 1 19961) . which calculate the FWHM (full-width 
half-maximum). Twenty percent of the peak intensity is related 
with the FWHM as W 20 = 1.524 FWHM. Finally, we corrected 
for the instrumental width using the expression 

AA = JaA 2 , -AA 2 . (2) 

y obs ms v ' 

The rotational velocities obtained from the widths of differ- 
ent optical lines are compared in Fig. 3. Although all veloc- 
ities are remarkably similar, a larger dispersion in the H^ vs 
[OIII]/15007A can be seen. The reason is that H^ has low signal- 
to-noise in most spectra and is heavily affected by stellar absorp- 
tion, so this line has not been used in the study. The fit is quite 
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Fig. 4. Comparison between local TFR obtained from DEEP2 data (0.1<z<0.3) in B-band, and I Verheijenl (1200 ll) local TFRs. The 
solid line represents the linear fit to the DEEP2 data (See text for definitions of the subsamples). 



good in the last panel, where [OIII]/15007A is compared ver- 
sus [OII]/l/L3727A, although only four objects have both lines 
in common. It is important to point out that a fifth galaxy with 
[OII]/Ll3727A velocity, 50 km/s higher than that obtained using 
[OIII]/15007A, was discarded because its spectrum shows that 
the object was an AGN or another active object. This leaves a 
final sample of 344 galaxies for which rotational velocities were 
calculated. 

3.2. Absolute magnitudes 

The absolute magnitudes were obtained after correcting for 
Galactic reddening, K-correction, and extinction. DEEP2 mag- 
nitudes are provided already corrected for G alactic reddening 
based on Schlegel, Finkbeiner & Davisl d!998l) dust maps. 

The K-correc t ion w as performed by adopting the 
Blan ton & Roweisl (120071) procedure (code version v4_l_4), 
which provides the K-correction from the B, R, and I-band 
photometry. The reliability of the K-corrections thus obtained 
depends mainly on the errors of DEEP2 photometric catalogue, 
which are lower than 0.2 mag. After K-correction, absolute 
magnitudes and its inverse variance were calculated from the 
luminosity distance corresponding to the measured redshifts, 
using the concordance cosmology. 

Finally, the absolute magnitudes were corrected for intrinsic 
extinction. This correction is basically dependent on inclination. 
Several works h ave studied this effect , providin g different 
correc tions. The iTullv & Fouqud d 19851) and the iTullv et all 
(1998) methods are the most common corrections used in TFR 
studies (e.g. lBohm eTaIll2004tlChiu. Bamford & Bunkerj|2007l) . 
In this work we have adopted the Tullv et al.l d 19981) method so 
as to compare it with the most recent TFR studies. According to 
this method, the extinction, A^, as a function of inclination, i, in 
the A - band, is defined to be: 



4"° = 7a log (a/b) 



A 1 ' 
with 

y B = -0.35(15.31 + M B ) 
y R = -0.24(15.91 + M R ) 
ji = -0.20(16.61 + Mi) 

where a/b is the galaxy major-to-minor axis ratio. 
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Fig. 5. Comparison between local TFR obtained from DEEP2 
data (0.1<z<0.3) in R-band (up) and I-band (down), with the 
I Verheiienl (l200li) local TFR of the 'RC/FD sample (without 
NGC 3992)'. The latter relation is determined from V max , and 
the galaxies included have V max > Vfi at . The solid line repre- 
sen ts the linear fit t o the DEEP2 data, while the dashed line is 
the IVerheii enl d200 lb relation in the respective band. 



4. Results 

4.1. Local Tully-Fisher relation 

The local TFR in the opt ical bands have been studied in severa l 
works over the years (e.g. lTully & Piercei2000HVerheijenl200lh . 
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Table 1. Best approximation between this work and lVerheiienl d2001l) parameters obtained in the sample (iv), by fitting V, 







This work 








Verheiien (2001) 


Band 


a 


h 


@~ total 


A 


B 




A 


B 


B 


0.160 + 0.062 


-0.113 ±0.003 


0.149 


1.42 ± 0.55 


-8.85 ± 


0.24 


1.24 ±0.82 


-8.5 ± 0.4 


R 


0.183 ±0.057 


-0.108 ±0.003 


0.144 


1.69 ±0.53 


-9.26 ± 


0.26 


1.54 ±0.82 


-8.9 ±0.4 


I 


0.219 ±0.055 


-0.105 ±0.003 


0.142 


2.09 ± 0.53 


-9.52 ± 


0.27 


3.05 ± 0.83 


-9.6 ± 0.4 



The slope and zero-point of the relation are very similar in all of 
them; however, a small difference in one of these parameters can 
be definitive for finding evolution. Therefore, a local TFR from 
the same data and in the same conditions as the high-z study 
must be established, for a reliable comparison. 

The lVerheiienl (120011) local relations are widely used in most 
recent evolution studies o f the TFR. For th is reason, our local 
fits are compared with the Verheiienl (12001 1) relations as a refer- 
ence. In this paper, different relations are derived as a function of 
the rotation velocity. The first one uses the parameter W' R r de- 
fined in iTullv & Fouqud d 19851) . as the rotation motion param- 
eter derived from the observed HI line width at 20% of peak 
intensity, corrected for random motions, transition from horned 
to Gaussian profile, and inclination. With optical lines, only the 
inclination correction is required because we assume Gaussian 
profiles, and, as stated above, random motions are less impor- 
tant in the optical. 

The other two relations are calculated from the rotation 
curve. One of them is derived from the maximum rotation ve- 
locity, and the other one from the velocity in the flat region of 
the rotation curve. Since these two vel ocities might not be the 
same, different samples are used in the VerheiienT d2001l) paper. 
In the present work, the velocity is calculated from line widths 
so that it does not correspond exactly with any of them. For this 
reason, the results of the local TF R obtained from DEEP2 are 
compared with all IVerheiienl (12001 1) relations derived from dif- 
ferent samples: 

(i) The 'RC/FDR sample' corresponds to all galaxies with 
resolved rotation curve, including the objects for which the ro- 
tation curve is truncated before reaching the flat region (only 

Vmax)- 

(ii) In the 'RC/FD sample', the galaxies whose rotation curve 
is truncated before flat region were eliminated, so it is possible 
to calculate V max and Vfl at . 

(iii) In the 'RC/F sample', only galaxies with V max =Vfl at are 
included. 

(iv) In the 'RC/FD sample (without NGC 3992)', this galaxy 
has been excluded from the 'RC/FD sample' because it was re- 
sponsible for a larger scattering in the relation. 

Our local TFR in the B, R, and I-bands was ob- 
tained from least-squares fitting to the DEEP2 data points 
in the redshift range 0.1<z<0.3. Since no significant change 
in the stellar populations of field galaxies was observed 
in this redshift range, we can assume that their TFRs 
are represent ative of the local ones. We followed other 
authors (e.g. iBamford. Aragon-S alamanca & Milva ng^-Jensenl 
2006; iNakamura et al.ll2006l:lvogtll2007l) by adopting loz(2V m „^ 
as the dependent variable in the fit: 



log(2V ffl „) = a + b M Q 



(7) 



where V max is the maximum rotation velocity, a the zero-point 
of the relation, b the slope, and Mq the absolute magnitude in the 
Q-band. This is the so-called inverse TFR, which is less sensi- 
tive to luminosity incompleteness bias ( Willick 1994; Schechter 
1980). 



The best fit to the data was found by minimizing the chi- 
squared statistics, 



X" 



2^ Wi (log(2V maxJ ) -a-b M QJ ) 2 



(8) 



where w ; = 1 /erf is the weight applied to each point according 
to the expression: 



_2 2 , r.2 _2 

°1 = °7og(2V„,„,,) + b °lf Q , t 



(9) 



Here, a-\ og (2v„ m ,i) and <tm Qi are the errors derived from the re- 
duction process. These weights are recalculated in each iteration 
with the new slope, b. 

The total scatter in the Tully-Fisher fit is calculated using 



total 



2^ Wi (log(2V mfl . M ) -a-b M QJ f 

i 



(10) 



In the Verheiienl (1200 ll) tables, their slope B and zero-point A 
correspond to the expression: 



Mq = A + B log(V) , 



(11) 



where log V is either log W' RI , log (2V max ) or log (2V/ to ), de- 
pending on the case. 

Figure 4 shows the c omparison of our local DEEP2 TFR in 
the B-band, with the four ! Verheiienl d200 lb samples. The smaller 
difference between both relations is obtained using the sample 
(ii) from Vfi at . However, when considering the three bands all 
together, the smallest difference is obtained using the sample (iv) 
by fitting V max . All th e param eters obtained are within the errors 
derived by Verheije^ d2001l) (see Table 1). Figure 5 shows the 
comparison between this sample and the DEEP2 TFRs in R and 
I-bands; therefore, the TFRs obtained from DEEP2 data in the 
redshift range 0.1<z<0.3 were adopted as our local TFR, since 
they were derived in a similar way to the higher redshift galaxies. 

4.2. Tully-Fisher relations in B, R, and I-bands 

With all corrections applied, the TFRs in B, R, and I-bands are 
obtained, dividing the redshift range 0.1<z<1.3 into 6 bins, as 
shown in Fig. 6 for the B-band. At the highest redshift, only the 
brightest galaxies can be seen, and since the scatter spreads the 
data points, the fit was performed by setting a constant slope, to 
study the evolution of the zero-point. In other words, a change in 
the slope cannot be distinguished from a change in zero-point. 
The fits were made in the same way as the local fit, but in this 
case no iterations were needed because the slope was considered 
invariant. 

As can be seen in Fig. 6, the zero-point of the TFR changes 
with redshift, since the difference with the local zero-point is 
larger as the redshift increases. For the last bin, corresponding 
to l.l<z<1.3, the difference is more than 2cr. It is important to 
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Fig. 6. Tully-Fisher relation in the B-band for 6 redshift ranges from 0.1 to 1.3. Thin and dashed lines represent, respectively, the 
local TFR and its 2<x uncertainty in the offset. The thick line is the weighted least-squares fit to the data, and the three dot-dashed 
lines are its 2cr uncertainty in the offset. 



note that the same zero-point change can be seen in both R and 
I-bands. A similar result for the lower redshift bins was found by 
IVogtl (120071) . where all fits had lower zero-points than the local 
TFR. Nevertheless, in her work, all fits were within 2<x, being 
0.9<z<l.l the last redshift range. She concluded that no TFR 
evolution was detected and that the change of the zero-point was 
instead attributed to objects with high equivalent widths. In Fig. 
6, our result is the same up to z - 1.1, but in the next redshift 
range (1.1 <z< 1.3), not only is the zero-point larger than 2<x un- 
certainty in the offset, but most galaxies (~ 85%) have their con- 
tinuum brighter than their local counterparts. Then, it is unlikely 
that this zero-point change can be claimed as due to emitters of 
high equivalent widths. Nevertheless, and given the data scatter, 



a clear-cut statement on the evolution of the TFR up to z= 1 .3 re- 
quires higher redshift data, or data with lower associated errors. 
This type of work can be carried out using the new generation 
of 30 meter-class telescopes or with the deeper observations of 
DEEP3. 

To quantify the changes in the intercept of the TFR, it is help- 
ful to work with the magnitude residuals resulting from the re- 
lation, that is, the difference between the absolute magnitude of 
one galaxy and the magnitude from the local TFR (for its veloc- 
ity). In Fig. 7, these residuals versus redshift are represented in 
the B-band. A tendency for galaxies to be brighter in the past 
can be seen, although at the level of the relation scatter. In ad- 
dition, we represented the weighted average to the data in the 
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Fig. 7. Magnitude residuals of the DEEP2 sample as a function 
of redshift in the B-band. Only the galaxies with a velocity er- 
ror lower than 5% have been represented, so that the figure is 
visually clear. The black points are the weighted average in each 
redshift range represented in Fig. 6. 

same redshift ranges as in Fig. 6. No change in magnitude for z 
< 0.6 can be seen. Then, the galaxies become brigher as redshift 
increases. 

Although the data could be fit by a non-linear function, since 
data dispersion is quite high, it is difficult to determine whether 
or not a non-linear fitting can really provide a more realistic ap- 
proximation, so we have chosen a linear fit, which provides a 
qualitative idea of the evolution found. The least-squares fit to 
the data in each band, weighted by the errors in AMq and ob- 
tained through error propagation, yields the following relations: 



AM B = (-0.13 ±0.12) + (-1.00 ± 0.15) z 
AM R = (0.00 ±0.12) + (-1.00 ± 0.15) z 
AM/ = (0.04 ± 0.12) + (-1.01 ± 0.16) z. 



(12) 
(13) 
(14) 



However, this procedure is only valid if the TFR slope does not 
change with redshift, because the magnitude range is not the 
same at low as at high redshift. 

In the three bands, the results are very similar, with a change 
in magnitude at z=l, AM= -1.1 ±0.2 in B and I-bands, and 
AM= -1.0±0 .2 in R-band. This result is in agreement of 
that fo und byjBamford, Arag on-Salamanca & Milvang-Jensenl 
(2006) and Bohm et al. (2004), although these studies were car- 
ried out using only th e B-band. In a more recent work by 
iBohm & Zieglerl (120071) . a luminosity evolution of (-1.22 ± 
0.40)z was found for galaxies with redshift < z >= 0.5. Taking 
the errors into account, this result would agree with our work, 
too. However, these authors claim that this difference is found 
only for low-mass galaxies, whereas the distant high-mass spi- 
rals are compatible with the local Tully-Fisher relation. In our 
case, the h igh-mass galaxie s are b righter than local galaxies, in 
contrast to Bo hm & Zieg ler (2007) work. 

5. Discussion: Tully-Fisher evolution and its 
causes 

A change in the TFR with lookback time represents a variation 
in luminosity, velocity, or both. A velocity change can happen 
if the mass of the galaxy disc changes with redshift, for a fixed 
luminosity. The mass can increase by galaxy merger processes, 



but this would only be noticeable if the mass of the merger is at 
least ~15% of the parent galaxy, w hereas in this case the merger 
would destroy the galaxy disc. Portinari & S om mer-Lar senl 
(2007) discuss the TFR evolution as predicted by cosmologi- 
cal simulations of disc galaxy formation, including hydrodynam- 
ics, and all the relevant baryonic physics (star formation, chemi- 
cal and photometric evolution, metal-dependent cooling, energy 
feedback). These simulations predict a growth of the disc mass 
by infall of halo gas and ongoing star formation. This would 
increase the rotation velocity of the galaxy, but at the price of 
increasing galaxy luminosity as well. Then, in this model, the 
evolution would shift the position of the galaxies along the TFR. 
Then, a change in luminosity seems more likely than a velocity 
change. Within this scenario, assuming pure luminosity evolu- 
tion, a change in the slope of the TFR would mean a different 
evolution for high and low-mass galaxies, while a zero-point 
TFR change would indicate a uniform evolution. 



The lPortinari"& Sommer -Larsenl (|2007|) simulations show 
a B-band luminosity evolution AMb = -0.85 magnitudes at 
z=l, with a constant TFR slope. This result is less than the 
evolution found in this work, although in the same direction: 
galaxies would have been brighter in the past, for a fixed ve- 
locity. Their models predict a luminosity fading and redden- 
ing caused by the ageing in the stellar populations, combined 
with the increase of the stellar mass and rotation velocity as 
stated above. One point in favour of this theory is the increase 
in the star formation rate (SFR) density from z= up to z=l-2 
observed in numerous works (see, for example, Guzman et al.l 
119971: iMadau. Pozzetti & Dickinsod[l998l iFeulner et al.ll2005l) . 
Although the mass locked in stars would still actually be grow- 
ing, their growth would have been greater in the past, affected 
by stellar populations ageing. Were this explanation correct, the 
effect should be more noticeable in blue-bands. However, in this 
work, the same change in luminosity in B, R, and I-bands has 
been found. 

The local TFR is colour dependent, in other words, the 
galaxies that rotate faster are redder. In Fig. 8 we repre- 
sent the colour R-I of our galaxies versus the rotation ve- 
locity and the local difference between R and I-band TFRs. 
For the three redshift ranges, all galaxies follow the colour of 
the local relation, although with larger scatter for higher red- 
shift galaxies. Then, the colour of a galaxy for a fixed veloc- 
ity does not change significantly with redshift. Also, were the 
slope of the TFR constant with redshift as other works sug - 
gest dBamford. Aragon-Salamanca & Milvang-Jensenl |2006), 
the colour of the Tully-Fisher relation would not change with 
the redshift. 



Finally. IWestera et al.l d2002l) studied the colour evolution of 
disc galaxies in two models of galaxy formation: the accretion 
model, and the collapse model. In both models, the intrinsic in- 
tegrated colour evolution R-I would be lower than 0.2 in the 
redshift range 0<z<1.2, within the dispersion shown in Fig 8. 
However, the change in luminosity found in our work can only 
be reproduced by the collapse model. The greatest colour evo- 
lution difference between both models can be found in V-K. In 
this case, while the collapse model presents a change in colour 
lower than 0.1, the colour evolution in the accretion model at 
z=1.2 would be 0.6. Then, the study of optical-infrared colours 
in TFRs can provide insight into the mechanism of formation of 
disc galaxies. 
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Fig. 8. R-I absolute magnitude versus rotation velocity in the redshift range 0.1<z<1.3. Solid line shows the difference between local 
R and I band TFRs. The crosses indicate the galaxies in the redshift range 0.1<z<0.3. The squares are galaxies with 0.3<z<0.8, and 
the diamonds represent the objects with 0.8<z<1.3. 



6. Conclusions 

This study investigates the Tully-Fisher relation and its evolu- 
tion with redshift in B, R, and I-bands. Using HST imaging 
data and DEEP2 spectroscopy taken with the Keck II telescope, 
galaxy characteristic parameters, such as rotation velocity and 
absolute magnitude, have been derived. The rotation velocity of 
a sample of 344 spiral galaxies in the Groth field were measured 
from optical lines widths. The objects span the redshift range 
0.1<z<1.3. In this work we demonstrated that all optical emis- 
sion lines can be used for determining disc-rotation velocities, 
if the signal-to-noise is high enough. Therefore, it is possible to 
extend the study of the TFR in the optical domain up to z=1.6. 

In addition, we used three different techniques for morpho- 
logical classification: comparison of the spectra with an ellip- 
tical/SO template, visual classification, and GIM2D. The best 
method was the visual classification, and we found an error of 
~25% in the selection, based on the Sersic index. 

Local Tully-Fisher relations were constructed from DEEP2 
data in th e reds hift range 0.1<z<0.3, and were compared with 
IVerheiienl d2001l) relations. Despite the different ways of deter- 
mining rotation velocity, the local relations derived in this work 
are consistent with his results, so our local TFRs have been 
adopted as low-redshift comparison patterns. 

The B, R, and I-band TFRs were derived for 6 redshift 
ranges. A least-squares fit of the data in each range was ob- 
tained, maintaining the slope constant, for studying the evolution 
of the zeropoints. We find prima facie evidence of evolution of 
the TFR, with a tendency for galaxies to be brighter in the past. 
The same zero-point changes can be seen in R and I-bands. We 
suggest that this evolution is due to an ageing of the stellar pop- 
ulations as consequence of a star formation decrease down to 
z=0. 

The colour and luminosity evolutions found in our work sup- 
port the collapse model rather than the accretion model of disc 
galaxy formation, showing that V-K colour evolution of TFRs 
can provide useful tools for studying galaxy formation. 
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